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ABSTRACT 

The yellow supergiant content of nearby galaxies can provide a critical test 
of stellar evolution theory, bridging the gap between the hot, massive stars and 
the cool red supergiants. But, this region of the color- magnitude diagram is 
dominated by foreground contamination, requiring membership to somehow be 
determined. Fortunately, the large negative systemic velocity of M31, coupled to 
its high rotation rate, provides the means for separating the contaminating fore- 
ground dwarfs from the bona fide yellow supergiants within M31. We obtained 
radial velocities of ~2900 individual targets within the correct color-magnitude 
range corresponding to masses of 12M & and higher. A comparison of these veloc- 
ities to those expected from M31's rotation curve reveals 54 rank 1 (near certain) 
and 66 rank 2 (probable) yellow supergiant members, indicating a foreground con- 
tamination > 96%. We expect some modest contamination from Milky Way halo 
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giants among the remainder, particularly for the rank 2 candidates, and indeed 
follow-up spectroscopy of a small sample eliminates 4 rank 2's while confirming 
5 others. We find excellent agreement between the location of yellow supergiants 
in the H-R diagram and that predicted by the latest Geneva evolutionary tracks 
which include rotation. However, the relative number of yellow supergiants seen 
as a function of mass varies from that predicted by the models by a factor of 
>10, in the sense that more high mass yellow supergiants are predicted than are 
actually observed. Comparing the total number (16) of > 2OM yellow super- 
giants with the estimated number (24,800) of unevolved O stars indicates that 
the duration of the yellow supergiant phase is ~3000 years. This is consistent 
with what the 12M and 15M evolutionary tracks predict, but disagrees with 
the 20,000-80,000 year time scales predicted by the models for higher masses. 

Subject headings: supergiants — stars: evolution — galaxies: stellar content - 
galaxies: individual (M31) 

1. Introduction 

The color magnitude diagrams (CMDs) of nearby galaxies reveal the details of stellar 
evolutionary processes, if our eyesight is keen enough. In recent years both the data and 
theory have made considerable advances. On the one hand, large format CCD cameras have 
allowed comprehensive photometry of the resolved stellar content of nearby galaxies, such as 
that of the Local Group Galaxy Survey (LGGS), which imaged those Local Group galaxies 
with active star formation (Massey et al. 2006, 2007a, 2007b). On the other hand, recent 
advances in stellar evolutionary theory have demonstrated the important role that rotation 
plays in the evolution of massive stars (see, for example, Maeder & Meynet 2000; Meynet & 
Maeder 2003, 2005). 

Consider the optical CMD of the Andromeda Galaxy (M31), shown in Figure [T] up- 
per. The unevolved stars are found on the left, in the section labeled "Blue Supergiants". 
However, this "blue plume" in such diagrams actually contain a mix of both unevolved main- 
sequence and more evolved blue supergiants (Freedman 1988); it also contains a smattering 
of Wolf-Rayet stars, the evolved descendants of the most massive O-type stars. Massey et 
al. (1995) have emphasized that the optically brightest stars in this region of the CMD are 
not, in fact, the most bolometrically luminous or massive — rather, the brightest stars are 
dominated by B- and A-type supergiants, while more massive (and luminous) O-type stars 
are fainter as most of their radiation lies in the far-UV. 
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The central portion of the CMD contains the yellow supergiants, and at the extreme 
right, the red supergiants. But, redwards of the blue plume caution is advised in our in- 
terpretation, as foreground contamination may dominate. We demonstrate this in Figure [T] 
lower by using the Besancon model (Robin et al. 2003) of the Milky Way to construct a 
theoretical CMD of the expected foreground contamination, using the same galactic coordi- 
nates as M31, and covering the same area. For the bright stars (V < 20) redwards of the 
blue plume (B — V > 0.4) the foreground contamination is > 70%. Note that the features 
in the M31 "yellow plume" area are well reproduced by the Besancon model. In Figure [5] we 
break down this foreground contamination into its various components. The disk popula- 
tions clearly dominates, but even among the brighter stars there will be some contamination 
by halo giants and sub-giants. 

Our group is engaged in the long-term process of characterizing the massive star popu- 
lations of nearby galaxies from one side of the CMD to the other. In order to relate this to 
evolutionary theory, we must succeed in two things: first, to be able to clean foreground stars 
from the sample, and secondly to provide a transformation of observed properties to physical 
properties. We have recently undertaken this for the red supergiants of M31 (Massey et al. 
2009); here we turn our attention to the yellow supergiants. 

1.1. Yellow Supergiants as a Magnifying Glass 

Yellow supergiants (F- and G-type) are extremely rare, as they represent a very short- 
lived phase in the evolution of massive stars. Their numbers and location in the H-R diagram 
(HRD) are very sensitive to the uncertain values of the mass-loss rates for massive stars, 
and how convection and other mixing processes are treated (Maeder & Meynet 2000). As 
Kippenhahn & Weigert (1990) put it, "[The yellow supergiant] phase is a sort of magnifying 
glass, revealing relentlessly the faults of calculations of earlier phases." 

Exactly how sensitive our expectations should be to the details of the models is illus- 
trated by comparing the various Geneva evolutionary models (Schaller et al. 1992; Charbon- 
nel et al. 1993; Schaerer et al. 1993; Maeder & Meynet 2001; Meynet & Maeder 2003, 2005) 
shown in Figure [3] and the corresponding lifetimes for the yellow supergiant stage given in 
Table [fl The metallicities shown span a range of 10, from sub-solar (z = 0.004, typical of 



1 Note that the 20M Q and 25M Q z=0.040 models of Meynet & Maeder (2005) were computed using 
a numerical simplification which resulted in the tracks turning back to the blue at too high an effective 
temperature; here we use the recomputed versions mentioned in Massey et al. (2009). We also include here 
newly completed z = 0.040 models for 12Mq and 15M® which include the effects of rotation. 
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the SMC) to solar (z = 0.020) to super-solar (z = 0.040, typical of M31); see Massey (2002) 
and references therein. Solid curves denote the latest models computed with an assumed 
initial rotation of 300 km s -1 ,the dashed curves are the latest models that include no initial 
rotation, and the dotted curves are the older, non-rotating models. Since different versions 
(rotating, newer non-rotating, older non-rotating) of the same mass track differ in luminos- 
ity, we have color coded them for clarity; the corresponding initial mass is shown in the same 
color. The two solid vertical lines denote the yellow supergiant region, which we define as 
having effective temperatures between 4800 K and 7500 K. 

We see that in most cases the models predict a short pass through the yellow supergiant 
region as stars evolve from the OB main-sequence stage (off the plot to the left) over to the 
red supergiant (RSG) region on the right. However, in some cases, such as the evolutionary 
track that includes rotation (solid curve) for 25M at Galactic metallicity (z = 0.020) the 
star then evolves back to the blue side of the HRD. At higher masses (40-60M Q ) the models 
predict that a star will double back to the blue side while a yellow supergiant. The lifetimes 
for the yellow supergiant stage given in Table [1] are all very short, typically a few tens of 
thousands of years. The main sequence lifetimes of these stars are a few million years, so, in 
general, the yellow supergiant lifetimes are on the order of 1% or less. In the extreme cases 
(lifetimes of the yellow supergiant phase of ~ 3000 yr for a 12M Q star with 15 Myr lifetime) 
it is of order 0.02%. We do see from Table [U however, exactly how sensitive these yellow 
supergiant lifetimes are to the exact treatments of the models. At Galactic metallicities, 
the models computed with an initial rotational speed of 300 km s _1 ("S3") indicate a 10 x 
shorter yellow supergiant phase than models with the same assumptions but with no initial 
rotation ("SO"). 

There are also clearly differences in the expected locations of the yellow supergiants 
in the HRD, particularly those of the highest luminosities. According to the models, for 
instance, at log L/L & ~ 5.8 dex we would expect to see some yellow supergiants at Galactic 
(z = 0.020) metallicities, but only the hotter ones (i.e., T e g >6000 K), while at lower 
luminosities we should find stars populating throughout the region. The duration of the 
yellow supergiant phase is not significantly shorter for a higher mass (luminosity) star than 
for a lower mass (luminosity) star, according to Tabled] — if anything the opposite is true — so 
simply determining the upper luminosity limit to yellow supergiants at various metallicities 
would be of great interest for comparison with the models, as well as seeing if the number 
of higher mass yellow supergiants really comparable to that of lower masses. 

We emphasize that such a test is new, and avoids some of the selection biases that may 
dominate comparisons of one population of massive stars to another, such as comparing the 
number of yellow and red supergiants, or the number of red supergiants and Wolf-Rayet 
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stars. Although those tests are invaluable, they require a thorough understanding of the 
completeness of the surveys of such objects. We are making progress towards gathering the 
information that make such further tests possible (see summary given by Massey 2009). 

1.2. Separating the Wheat from the Chaff 

We have previously found that it is straightforward to separate (extragalactic) red su- 
pergiants from (foreground) red dwarfs and giants by using a two-color diagram (Massey 
1998a, Massey et al. 2009), but this discrimination does not extend to the yellow super- 
giants. In Figure H] (upper) we show the intrinsic color lines from FitzGerald (1970). The 
dwarf sequence is shown in green, and the supergiant sequence in red. We have reddened the 
supergiant sequence slightly, by E(B — V) = 0.13, corresponding to the typical reddening of 
a massive star in M31 (Massey et al. 2007a); this slight adjustment down and to the right is 
shown by the short bar in lower left. Superimposed on this two-color diagram are the stars 
brighter than V = 18.5 from the Local Group Galaxies Survey (LGGS) photometry given 
by Massey et al. (2006) for stars in M31. 

For the mid F-type stars, there is a separation in the sense that for a given B — V a 
supergiant will have a larger U — B. For early A- type supergiants this trend is reversed; i.e., 
supergiants will have a more negative U — B. However, as is clear from this plot there is 
little or no separation in the two-color diagram for supergiants of late A or early F, or for 
supergiants of late F through early K. Examination of other colors (using model atmospheres) 
failed to identify any more suitable diagnostic two-color tool. 

It is also clear from distribution of points in this figure that the majority of blue stars are 
expected to be supergiants, while the majority of stars of later types may well be foreground, 
as we argued above. We illustrate this further in Figure H] (lower), where we have used the 
Besancon simulation of the Milky Way (Robin et al. 2003) to show the expected distribution 
of foreground stars in such a diagram. The comparison between the two sides of the figures 
reveals where we expect foreground stars to dominate. 

Given the lack of two-color discriminants, what then are our options to identify a rel- 
atively complete sample of yellow supergiants that can be used to test the stellar evolution 
models? Uncertain distances and reddening complicate the analysis of a galactic sample. In 
her compilation of 949 supergiants known in the Milky Way, Humphreys (1978) lists just 
21 supergiants of spectral type F0-G9 I (i.e., about 2%). Even so, several of these have 
dubious cluster memberships and, hence, uncertain luminosities. In the Magellanic Clouds 
there are less than a dozen spectroscopically confirmed F and G supergiants (Oestreicher & 
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Schmidt-Kaler 1999), although this deficiency may be largely due to the lack of adequate 
spectroscopic studies in the correct magnitude /color range, a situation we ourselves hope to 
remedy in the not too distant future. 

We have concluded that M31 provides the best laboratory for conducting such tests 
at present, as its kinematics allow us to overcome the problems posed by foreground con- 
tamination of its CMD. M31 possesses a large negative systemic velocity (~ —300 km s _1 ) 
and a high rotation rate (240 km s _1 ), making it relatively straightforward to demonstrate 
membership based on radial velocities. Gilbert et al. (2006) and Koch et al. (2008) similarly 
used radial velocities to separate M31's red giant members from foreground contamination. 
The recent LGGS photometry of M31 (Massey et al. 2006) provides the means for selecting 
candidate stars for radial velocities, and for transforming intrisic colors to physical properties. 

In this paper we conduct a census of yellow supergiants in M31, establishing membership, 
determining physical properties, and making comparisons with the current generation of 
evolutionary tracks. In § [2] we describe our data and reduction. In § [3] we illustrate the 
process by which we separated foreground dwarfs from M31's yellow supergiants, and in §H] 
we present the comparison of our results to the current evolutionary tracks. We provide 
discussion and lay out our thoughts for future work in § [51 

2. Observations and Reductions 

In order to separate yellow foreground dwarfs from the yellow supergiants, we used the 
Hectospec 300 fiber spectrograph on the 6.5-m MMT telescope to obtain radial velocities 
for ~ 2900 stars. In this section we describe the sample selection, data acquisition, and 
reductions. 



2.1. Sample Selection 

Our sample of objects was selected to have V < 18.5 (roughly corresponding to log L/L Q ~ 
4.4) to provide adequate signal-to-noise for our spectroscopy. The color range was originally 
restricted to U — B > —0.4 with 0.4 < B — V < 1.4 in accordance with the range of 
B — V for which dwarfs and supergiants cannot be distinguished photometrically; however, 
we found that relaxing the color range to 0.0 < B — V < 1.4 added only a small percentage 
of additional stars. 

Our other selection criteria was based upon the need to be able to distinguish the M31 
yellow supergiants from foreground disk dwarfs based upon their radial velocities. Inspection 
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of the atlas of Galactic neutral hydrogen by Hartmann & Burton (1997) shows that M31 (at 
I = 121.2°, b = —21.6°) clearly stands out from the Galactic clutter at a radial velocity of 
-150 km s _1 (see Hartmann & Burton 1997, p. 87), but is somewhat confused by -100 km s _1 
(pp. 95, 97). We therefore decided to restrict our observations to those areas in M31 where 
the radial velocities should be < -150 km -1 . For computing the expected radial velocity 
corresponding to a position in M31, we used Rubin & Ford (1970), the seminal paper on the 
rotation of M31. A least-squares linear fit to the radial velocities of the HII regions yields 
the expected radial velocities Vel ex pcct 

Vel expcct = -295 + 241.5(X/i?) 

where X is the distance along the semi-major axis, and R is the radial distance of the object 
within the plane of M31. We found that this approximation works well, producing good 
agreement with the more complex two dimension velocity field (Sofue & Kato 1981), and 
with other recent approximations (Hurley-Keller et al. 2004). The radial velocities of red 
supergiants in M31 also agrees well with this simple relationship (Massey et al. 2009). 

The result of this selection criterion is that not all of the area surveyed in the LGGS 
was included in our sample: stars along the south-west half of the semi-major axis of M31 
(X/R ~ —1) will have radial velocities of ~ —550 km s _1 , but stars along the north-east half 
of the semi-major axis (X/R ~ 1) will have radial velocities of —50 km s _1 , more positive 
than our selection criterion of —150 km s -1 . The distribution of stars in our sample is shown 
in Figure O where the "jaws" in the north-east (upper left) are due to stars near the semi- 
major axis having smaller R values then those seen along the edge of the disk. In all our 
sample covered 1.6 deg 2 of the 2.2 deg 2 of the LGGS survey (i.e., about 73%). 

We were concerned that a few legitimate F or G supergiants might be too bright to 
be included in the LGGS photometry, particularly due to saturation in the R filter around 
R ~ 15.5. We therefore supplemented the LGGS data with 163 bright (V < 16.0) stars from 
the survey of Magnier et al. (1992). In order prevent mixing bright and faint stars in the 
same observations, we then divided our overall catalogs into a "bright" catalog (349 stars 
with V < 16.0) and a "faint" catalog (3994 stars with 15.5 < V < 18.5), with 61 objects 
in common. Due to the constraints of the fiber configurations, not all the objects could be 
observed, but we did manage to observe 68% of the 4282 catalog targets. 

2.2. Spectroscopy 

Hectospec is a 300 optical fiber fed spectrograph (Fabricant et al. 2005) on the MMT 
6.5-m telescope. Observations are obtained in an innovative queue mode where the observers 
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are the astronomers who were awarded time, but with the observing program for a given night 
determined by a queue manager. The observers are ably assisted by one of two professional 
instrument operators, in addition to the telescope operator. This "collective" approach 
spreads out the effects of poor weather throughout the observing season, reducing the impact 
on any one program. Our observations were all carried out on eight nights during October 
2007 and one night in November 2007. The 600 line mm -1 grating was used, providing a 
dispersion of 0.55 A pixel -1 and a (5 pixel) spectral resolution of 2.8 A. The wavelength 
coverage extended from 4550-7050 A. 

The fiber configuration files were constructed prior to the observations. The instrument 
has a 1° field of view, a reasonable match to the 3-4° angular extent of the optical disk of M31 
(Hodge 1981) . Our observations consisted of observing a single configuration for each of the 
five fields containing the brighter stars ("Brt" fields), and multiple configurations of three 
fields containing the fainter stars ("Fnt" fields), as listed in Table [2j Observations each "Brt" 
configuration consisted of 3 consecutive exposures of 10 minutes each, while observations of 
the "Fnt" configurations consisted of 3 consecutive exposures of 15 minutes each. In the 
end, we obtained 3116 spectra of 2901 of our catalog targets. 

Owing to the logistics of the queue observations, calibration exposures (flat field and 
He-Ne-Ar) were taken only in the afternoon, and subsequent to that the grating might have 
been tilted to a different angle to accommodate other programs, or even removed entirely 
and then replaced for our observations. We deal with this complication as described below. 

We also needed observations of stars that could serve as templates for the cross-correlation. 
Since the 600 line mm -1 is not commonly used, we obtained our own observations of three 
radial velocity standards, HD 196850, HD 194071, and HD 213014. HD 213014 was observed 
on four different nights. 

2.3. Data Reduction 

The data were all reduced using the "hectospec" IRAF@ package, designed specifically 
for this instrument (Mink et al. 2007). The data were bias-subtracted, trimmed, and a bad 
pixel extrapolation was performed using pre-existing bad pixel maps. The flat field exposures 
were extracted for each fiber and normalized in order to make the pixel-to-pixel corrections. 
The He-Ne-Ar arc exposures were then extracted, and used to make a dispersion solution. 



2 IRAF is distributed by NO AO, which is operated by AURA under cooperative agreement with the NSF. 
We appreciate the on-going support of IRAF by the volunteers at the IRAF help "desk" , [http://www.iraf.netl 
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A sixth order chebyshev function was used for this, resulting in typically 0.04 A residuals. 

The program exposures were extracted, using the dome flat field exposures as reference, 
and wavelength corrected. As the He-Ne-Ar exposures were obtained in the afternoon, and 
the grating might even be removed and reinserted before the program exposures, a zero-point 
shift in wavelength was determined for each of the M31 spectra using the O I A5577 night sky 
line. (As explained below, no such correction could be made for the very short exposures of 
the radial velocity standards.) Consecutive exposures of each M31 configuration were then 
summed, after cosmic rays were first identified and removed by comparing a median of the 
exposures to the individual exposures. 

For sky subtraction, each fiber had to be first be corrected for its own wavelength de- 
pendent throughput, using either the dome flat exposures, or, preferably, twilight exposures, 
if the latter had been obtained. Each M31 configuration contained both preselected "clean" 
sky positions plus random locations that might prove clean enough to be used as a measure 
of the sky. For each program spectrum 6 of these sky spectra were selected. These were 
taken from positions nearby on the array in order to reduce any scattered light component. 
These skys were then used to construct an average sky for subtraction using the Singular 
Value Decomposition method (Mink & Kutz 2001). No sky subtraction was needed for the 
bright radial velocity standards. 

3. Analysis 

3.1. Observed Radial Velocities 

The radial velocities of our 2901 program objects were obtained through a cross-correlation 
with suitable radial velocity standards. All cross-correlations were performed in the IRAF 
package "fxcor" , which computes radial velocities via Fourier cross-correlation using a Gaus- 
sian to find the center and width of the calibration peak, following the method of Tonry & 
Davis (1979). Before the cross-correlations were computed, each individual spectrum had 
its continuum removed by first normalizing and then subtracting 1.0. The normalization 
was done interactively, using a 12th order cubic spline fit to the continuum. The cross- 
correlations were restricted to the range 4750-5550 A and 5600-6800 A in order to avoid the 
strongest night-sky emission. 

As mentioned above, although the [O I] A5577 emission line was adequate to correct the 
program spectra for the wavelength zero-point (necessitated by the logistics of the queue) 
there was, ironically, no equivalent way to correct the radial velocity "standards" . In order 
to resolve this issue, we used two velocity templates, created for a related project, generated 
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from 270 line mm -1 Hectospec observations (Caldwell et al. 2009). The templates were 
created by first deriving an initial velocity from library templates (typically a K giant star) 
for spectra selected from a catalog of M31 star clusters. The spectra were shifted to zero 
velocity and sorted into crude spectral bins (F- and K-type) at which point the best spectra 
from each type were then combined to make new templates and the whole process repeated. 
When we cross-correlated our "rv- standards" with these templates, the velocities produced 
by the two templates varied by < 1.1 km s _1 for all six observations. A third template, 
based upon early-type spectra (A-type) gave a velocity ~5 km s _1 more positive, and we 
chose not to use it n. 




In order to obtain suitable velocities for the radial velocity standards the relative ve- 
locities of the standards to the templates were obtained via cross-correlation and corrected 
to heliocentric velocities. We were very pleased to find this process yielded velocities that 
agreed with the IAU adopted velocities to better than 5 km s _1 . This indicates that if we 
had simply adopted the standard velocities, with no corrections, the resulting error would 
not have been large enough to affect our results. 

As a final check of our newly adopted velocities, the standards were cross-correlated 
against each other. HD194071 and the four observations of HD213014 all produced veloc- 
ities within 0.2 km s _1 of the standard value value, and within 0.1 km s" 1 of each other. 
HD 196850, however, consistently produced velocities 0.5 km s _1 more positive. Although 
this value is not significant compared to the expected radial velocity difference between M31 
yellow supergiants and galactic yellow dwarfs, it still prompted us to reject our spectrum of 
HD 196850 as a radial velocity standard template. We then performed cross correlations of 
the five spectra of the radial velocity standards against each of the 3116 spectra of the 2901 
program objects. 

There were four objects for which our cross-correlations initially failed: J004101. 24+410434. 6, 
J004129.31+40502.9, J004203.63+405705.8, and J004459. 11+412732. 7. The spectra of all 
four objects were examined and strong nebular emission in the region of Ha was evident. 
When this region was excluded good correlations were found for J004101. 24+410434. 6, and 
J004129. 31+40502. 9, although their associated error was about twice as large as the average 
value of the other program objects. 

We show in Figure [6] examples of two of our spectra, and their resulting cross-correlation 



3 We did find that the low-resolution templates produced problems when measured directly against our 
M31 observations obtained at higher dispersion, giving inconsistent results on the one field repeated on 
two nights. The differences were strongly correlated with velocities, which we interpreted as being due to 
cross-correlating two spectra of very different dispersions and velocities. 
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functions and corresponding fits. These were selected to have V magnitudes typical of the 
median value in our entire sample, 17.1. 

We list in Table [3]the observed radial velocity of each of our 2889 objects, some of which 
had multiple observations. What are the corresponding errors? Each observed radial velocity 
vel bs is the average of cross-correlating the star's spectrum against that of the five radial 
velocity standard spectra. The standard deviation of the mean of these is always very small, 
of order a few tenths of a km s _1 . However, the signal-to- noise in the standard templates 
is extremely high, and the standards are of similar type, so the differences in the velocities 
obtained from cross-correlating the same spectrum against these different standard spectra 
will certainly underestimate the true uncertanty. Tonry & Davis (1979) instead introduce 
the r parameter, which is the ratio of the peak of the correlation function to its noise. We 
can estimate the relationship between r and the error by using measurements obtained for 
stars observed multiple times. The 128 objects in the Brt5-1 field were observed on two 
separate nights, and in addition, there were 28 objects in common between the "bright" and 
"faint" fields, plus 59 objects observed in multiple "faint" fields. We find that in general our 
velocity error (in km s _1 ) is given by 

Vel err = 2.3 + 11.5/(1 + r), 

where the functional form reflects both errors in the standard star velocities and in uncer- 
tainties due to the cross-correlation (see Tonry & Davis 1979). The typical (median) r value 
for our data is 33, corresponding to a 2.6 km s _1 error. 

As a further check, in Figure [7] we compare the velocities obtained from the two obser- 
vations of the Brt5-1 field. We find that there is an excellent match between the observations 
over all the velocities we've sampled. As mentioned above, this was not the case when we 
originally used the low- dispersion templates, where we saw a velocity-dependent problem. 
For all future purposes in the paper, the velocities produced during multiple observations 
were averaged to yield one Vel b s for each program object. It is clear that a few stars have 
variable velocities based on this comparison; we expect these stars to be binaries. If two 
observations of the same star differed by 10 km s^ 1 or more we note this fact in Table [31 

Since we will be using the supergiants we identify to test if the models predict the same 
relative number with luminosity, we should understand what magnitude biases exist, if any, 
in the subsample of targets for which we obtained radial velocities compared to the parent 
sample. In Figure [S] we show the distribution of magnitudes in the complete sample of 4282 
objects which met our original criteria (magnitude, color, and location) in black. The red 
(dashed) histogram shows the sample of 2989 objects for which we successfully obtained 
radial velocities. We show this on a log plot so that a linear difference corresponds to a 
percentage; i.e., a difference of 0.3 dex is a factor of 2 regardless of the absolute numbers. 
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We find very good agreement, with only a slight tendency towards having obtained radial 
velocities for a larger percentage of brighter stars than fainter ones. This will introduce 
a slight bias in our final results in that higher luminosity stars should be proportionally 
over- represented in our sample. 

3.2. Identifying the Supergiants 

Now armed with the observed radial velocities of our 2899 program objects, two main 
steps are necessary to identify the M31 yellow supergiants. First, foreground dwarf con- 
taminants must be eliminated by determining which objects' velocities match that expected 
from M31's rotation curve and, second, once a list of M31 members has been obtained, any 
non-stellar objects (i.e., small clusters) must be removed. 

In the notation of Rubin & Ford (1970), the radial velocity V r of an object in a disk 
galaxy can be approximated by V r = Vq + V(R) sin£ cos 8 where Vq is the systemic radial 
velocity, £ is the angle between the line-of-sight and the perpendicular to the plane of the 
galaxy, V(R) is the rotation velocition within the plane at a radial distance R, and cos# = 
X/R where X is the position along the major axis. This would be a linear relationship only 
if the rotation curve was absolutely flat i.e. V(R)= const. However, as mentioned in § 12.11 
this is a good approximation for M31. We used the same fit, as described there, to compute 
what the radial velocity of each of our program objects should have been (given its position 
on the sky) were it an M31 member. The difference between our V Q b s and V exp ect for all 
objects was then calculated and plotted against V expe ct, as shown in Figure [9j 

In Figure[9J all the objects whose velocities correspond with M31's rotation curve should 
lie along the zero point of the Y axis. The left hand side of the plot, where the expected 
velocities are highly negative, represents the SW portion of M31, rotating towards the Milky 
Way. Thus, the strong diagonal band represents the foreground dwarfs (objects with essen- 
tially zero radial velocity). As can be seen, the M31 members can easily be distinguished 
from the foreground dwarfs on the left side of the plot and the distinction becomes increas- 
ingly hazy as we move along M31's semi- major axis. We therefore distinguish between two 
"ranks" of M31 members. Rank 1 objects are those which we can say are "nearly certain" 
M31 members, whereas rank 2 objects we consider to be "probable" M31 members. All 
objects with an expected velocity < —280 and difference < 60 were labeled as rank 1, and 
those with expected velocity > —280 and difference < 60, expected velocity < —340 and 
difference < 180, or expected velocity < —440 and difference < 220 were labeled as rank 2. 
These distinctions are displayed in Figure [91 
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As a result of this classification we were left with a list of 56 nearly certain and 71 
probable M31 members. There still existed the possibility, however, that some of these M31 
members were not stars, but small clusters. This was addressed by measuring the size of 
the objects as compared to nearby stars. We did this using two methods, both utilizing the 
images of the LGGS. Initially the 127 potential supergiants were examined with the Source 
Extractor (Bertin & Arnouts 1996), using the V-band images, and the results were then 
compared to an independently conducted manual check on the R-band. The manual check 
was completed by measuring the FWHM of each potential supergiant and its neighbors in the 
IRAF package imexam. The agreement between the manual and Source Extractor method 
was complete with each yielding seven previously known (Galleti et al. 2007) clusters in our 
sample: Mag-237751, J004345.23+410608.5, J004446.42+412918.3, J004545.58+413942.4, 
J004356.46+412203.3, J004358.15+412438.8, and J004403.98+412618.7. 

Once the foreground contaminants and M31 clusters had been removed from our sample 
and our rank designations applied, we found that, out of our original sample of 2899 objects, 
we are left with 54 rank 1 and 66 rank 2 yellow supergiants. This corresponds to a foreground 
contamination between 96-98%! The velocity information for these 120 objects, as well as 
the seven M31 clusters (listed at the end) is summarized in Table SJ 

How clean is our remaining list of candidates? As discussed in § 11.21 we expect that our 
original sample to contain some small fraction of Milky Way halo stars in addition to the 
numerous foreground disk stars. Although the radial velocities are effective at eliminating 
the disk contaminants, they will be less effective at weeding out the halo stars. Of the 
~ —300 km s _1 systemic velocity of M31, about two-thirds of it is the reflex motion of the 
sun: equation 4 of Couteau & van den Bergh (1999) implies a reflex motion of —178 km s -1 . 
Stars from the halo are therefore likely to reflect this solar motion. If the halo's velocity 
dispersion is 130 km s _1 (Binney & Merrifield 1998), then 3a velocities could extend all the 
way to -570 km s~ l . If there was a significant number of such stars in our sample, then some 
contamination could occur. 

In Figure [10] we compare our observed velocities (dashed, red histogram) with that 
expected for foreground stars according to the Besancon model (solid, black histogram). 
We see that virtually all of the stars with observed velocities greater than —175 km s _1 
are foreground objects, as expected. At —200 km s _1 about half of the objects should be 
foreground, and half M31. More negative than this, the M31 population dominates, with 
increasingly little foreground contamination. The detailed output of the Besancon model 
shows 7 foreground stars with radial velocities < —300 km s _1 , and 1 foreground star with 
a radial velocity < —400 km . 

We can also estimate this contamination independently using the Bahcall & Soneira 
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(1980) model which predicts about 40 halo dwarfs per square degree within the magnitude 
and color range we use in the direction of M31. As for halo giants, the Bahcall & Soneira 
(1980) model overestimates their number (Morrison 1993), but Heather Morrision (private 
communication) has been kind enough to calculate that there should be about four such 
objects in a square degree seen towards M31 in this same magnitude and color range. The 
effective area in our survey is about 1 deg 2 , and in it we observed 68% of the stars in our 
original sample, so we might realistically expect about 30 halo stars in our spectroscopic 
sample. Assuming the radial velocity distribution is Gaussian, this leads to the identical 
results as above: there should be 5 halo stars with velocities more negative than —300 km 
s" 1 (0.94a) and 1 below -400 km s~ x (1.8a). 

So, simply noting that the rank 1 stars all have have velocities more negative than —280 
km s _1 , we expect at most 8 stars in this sample of 54, or 15%. However, this is likely an 
overestimate of the contamination, as it ignores the additional information gained by the 
star's position and therefore its expected radial velocity, V CX pcct, were it an M31 member. To 
be a rank 1 candidate, a star has to be in the right location within the M31 field such that its 
V bs corresponds to its V ex pect (which, although still possible, is improbable). We therefore 
refer to the rank 1 candidates as "nearly certain" . The situation for the rank 2 candidates 
is harder to evaluate, as a few of these have radial velocities as positive as —100 km s _1 , 
but their location in the Figure [9] leads us to consider them "probable" yellow supergiants. 
Contamination by foreground stars of half of this sample would not surprise us, however. 

Of our original 2899 program objects, four had previously appeared in the literature as 
spectroscopically confirmed members of M31. J003745. 26+395823. 6 appeared in Humphreys 
(1979) as IV-A207 and was classified as F5 la, J004101.24+410434.6 can be found by the 
name OB69-46 in Massey (1998a) as a red supergiant (based on earlier photometry), and 
J004129.31+405102.9 is listed as OB22A in Humphreys et al. (1990) and classified as a F8- 
G0 la. We categorized all three as rank 1 yellow supergiants. Additionally, Humphreys et 
al (1988) lists J004101. 55+403432.3 (as III-R61) as a M31 RSG (K5 I) candidate. However, 
we are forced to conclude that this star is actually a foreground dwarf, given its —41 km s _1 
radial velocity (at its position in M31, the expected radial velocity is —464 km s _1 . We have 
double checked our previous cross-identification of this star, and it matches the one shown 
on the Humphreys et al. (1988) finding chart; the colors B — V = 0.69 and U — B = 0.02 do 
not correspond to mid K-type, and the spectrum is clearly of earlier type, with the Balmer 
lines prominent, along with strong lines of Mg I and Fe I. 
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3.3. Physical Properties: Transforming the Photometry 

In order to use our list of M31 yellow supergiants to test the current stellar evolutionary 
tracks, it is necessary to determine their effective temperatures and bolometric luminosities. 
For this, we will transform each star's photometry, as described below. 

The B and V photometry for all but one of our stars came from the LGGS; that of the 
single "bright" star, Mag-253496, is taken from Magnier et al. (1992) and adjusted by the 
small correction found by Massey et al. (2006). We apply a constant reddening correction 
E(B — V) = 0.13 based on the median value found for early-type stars in M31 by Massey et 
al. (2007a), and is in accord with the color excess derived from spectral types of a handful of 
O-type stars (Massey et al. 1986). The reddening of individual stars can readily differ from 
this (by several tenths) but in this part of the CMD reddening-free indices such as Q are 
degenerate with T e ff, as one may notice in Figure HI the reddening vector is nearly parallel 
to the supergiant sequence. 

The problem now becomes how to best translate these dereddened colors into effective 
temperatures. Flower (1996) and Kovtyukh (2007), both present empirical effective tem- 
perature scales that include F- and G-type supergiants. The Flower (1996) data is drawn 
from the literature, while Kovtyukh (2007) performs his own analysis on spectra of Galactic 
supergiants. The problem is that both of these studies are based upon samples of Galactic 
supergiants, while we expect the metallicity of our M31 stars to be about 2x solar based 
upon HII region abundances (i.e., Zaritsky et al. 1994). We therefore have decided to instead 
use the Kurucz (1992) "Atlas9" model atmospheres to provide the transformations. 

In Figure [TT] we compare the Altas9 models for solar metallicity with the two empirical 
calibration^. For this comparison, we have used the lowest two surface gravities for each 
temperature model computed by Kurucz (1992), i.e., logg = 0.0 and 0.5 for T e fj < 6000 K, 
\ogg = 0.5 and 1.0 for 6250 K< T eff < 7500 K, and hgg = 1.0 and 1.5 for 7500 K< T eff < 
8500 K. For the purposes of comparing our data to the stellar evolutionary models, we will 
restrict ourselves only to our defined yellow supergiant effective temperature range: 7500 K 
to 4800 K (logTcfj = 3.875 and 3.681). This region is indicated by the two horizontal lines 
in Figure [TTJ We can see that over this temperature range there is substantial agreement 
between the models and the empirical calibrations. 

As mentioned above, one advantage of using the Atlas9 models is that we can fine- 
tune the transformations to an appropriate metallicity, although we will find below that this 



4 We have used the intrinsic colors of FitzGerald (1970) as a function of spectral types in order to assign 
(B — V)o values to the spectral types given by Kovtyukh (2007). 
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correction is very slight. In keeping with Zaritsky et al. (1994) we assume an M31 metallicity 
of 2x solar (see further discussion in Crockett et al. (2006) and Massey et al. 2009). Using 
the Atlas9 model with the most similar metallicity (1.6 x solar) we compute a relationship 
between (5 — V) and logT eff as follows: 

logTcfj = 3.913 - 0.3512(5 - V) + 0.2692(5 - V) 2 - 0.1108(5 - V)jj 

To keep the relationship accurate over the range of temperatures in which we are the 
most interested, we restricted the fit to models that just bracketed the temperature range 
above (i.e., Kurucz models with 4750 K< T c r < 8000 K), and therefore this relationship 
is only applicable for 0.03 < (5 — V)q < 1.26. Since we adopted an average value of 
5(5 — V) — 0.13, this then corresponds to 0.16 < B — V < 1.39, a good match to the 
0.0 < B — V < 1.4 of our sample. Only six stars of our 120 yellow supergiants have colors 
bluer than B — V = 0.16. Although we will not include these in our comparison to the 
evolutionary models in the next section, we nevertheless would like to include these in the 
H-R diagram, and so we will adopt the following approximate transformation for these stars: 
logTeff = 3.934 - 0.549(5 - V) . 

The bolometric corrections are relatively modest for yellow supergiants (a few tenths of 
a magnitude), and we derive the following relationship using the Atlas9 models; the results 
are valid btween 4750 K to 9500 K, and, hence, applicable to our complete sample: 

BC = -251.54+ 130.763 logT efr - 16.9934(logT efr ) 2 

How much difference does adopting 1.6 x solar metallicity make? The difference is slight: 
for a star with (5 — V)o = 0.6 (roughly corresponding to 6000 K) we would derive logT c ff = 
3.775 using the 1.6 x Atlas9 relationship given above, while we would derive logT e ff = 3.769 
from an analogous relationship derived from the 1.0 x models. The difference, 0.006 dex, is 
negligible. 

We give the derived physical properties in Table [5] for the 120 probable supergiants 
sorted by bolometric luminosities. A distance modulus of 24.40 (0.76 Mpc) was adopted, 
following the discussion in van den Bergh (2000). 

3.4. Membership Re-examined 

One of the critical and interesting properties we are attempting to determine is the 
upper luminosity limit for yellow supergiants, and we were struck by the fact that the ten 
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most luminous stars in Table [5] were all of rank 2. We argue above that we expect some 
(minimal) contamination by foreground (halo) giants in our sample, as these stars would have 
radial velocities characteristic of the reflex motion of the sun. In particular, the three most 
luminous possible supergiants stand out as extraordinarily bright. None of the three have 
extreme radial velocities (which is in part why they are rank "2"): Mag-253496 has a Vel b s 
of -227 km s" 1 , J004251. 90+413745.9 has a Vel obs of -210 km s" 1 , and J004618. 59+414410.9 
has a Vel b s of -150 km s 1 , as can be seen in Table HJ We were therefore very keen to 
confirm or refute their membership in M31. The wavelength range of our own spectra had 
been optimized for radial velocities and where Hectospec has good throughput, and did not 
include the various good luminosity indicators in the far blue (see below) nor the O IA7774 
triplet in the far red. 

The Kitt Peak Director was sympathetic to our plight and arranged follow-up spectra to 
be obtained for these three stars plus comparison spectral standards. The data were obtained 
as part of engineering time on the 3.5-m WIYN telescope with the Hydra fiber spectrometer 
in the far red, and on the Kitt Peak 4-m Mayall telescope with the RC spectrograph (in the 
blue). These spectra convinced us that all three of these are actually halo stars, and not M31 
supergiants, as we argue below; they also eliminated another rank 2 star from membership, 
and confirmed membership of other (both rank 1 and 2) stars. 

The WIYN Hydra spectra were obtained on 11 August 2008, and consisted of an 
1800 sec exposure and a 1300 sec exposure; the latter was ended due to clouds. The 
setup included Mag-253496 and J004251. 90+413745. 9 and a number of other stars, but not 
J004618. 59+414410. 9. The (2.8 pixel) resolution was 4.0 A. The wavelength range included 
the O IA7774 triplet, known to have a strong luminosity effect in F-type supergiants (Osmer 
1972) due to non-LTE effects, exacerbated by sphericity and the large mass outflows found 
in supergiants (Przybilla et al. 2000). We include in Table what we find for this line. The 
two most luminous supergiant candidates have little or no O I A7774, arguing they cannot be 
supergiants, and we identify another non-supergiant among the rank 2 objects we observed. 
Strong O I A7774 is found for several of the other rank 1 and rank 2 candidates. 

The 4-m RC Spectrograph spectra were obtained on 4 September 2008, and consisted 
of a 3x900 s exposure of Mag-253496, and 3x500 s exposures of J004251. 90+413745.9 and 
J004618. 59+414410.9. The wavelength region included 3880-4600 A at a (2.2 pixel) reso- 
lution of 1.6 A. Several spectral standards were also observed to provide guidance in inter- 
preting the data; these were supplemented by similar data obtained at a later data with the 
Kitt Peak 2.1-m Goldcam spectrometer at a similar dispersion. These data show that Mag- 
253496 is roughly of spectral type G8, based upon the Fe I A4143/ H5 and Fe I A4045/H5 
ratios (Keenan & McNeil 1976). J004251.90+413745.9 and J004618.59+414410.9 are of ear- 
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lier type, between F5 and G2, based upon the strength of the G-band compared to H7, and 
we adopt an F8 type. 

For the G8 star, Mag-253496, we find that the ratio of Sr II A4077 to the Fe II/III blend 
at A4063 is about 1, typical of an G8 III; in a supergiant, this ratio would be considerably 
larger (e.g., Keenan & McNeil 1976). Similarly, the ratio of the Fe, SrII blend at A4216 to the 
Ca II A4226 is quite small, consistent with that of a dwarf or a giant, but not a supergiant. 
Combined with the lack of O I A7774, we therefore conclude this star is a Milky Way halo 
giant, and not an M31 supergiant. For the F8 stars, we find that the strengths of Till 
A4444 and Mg II A4481 to be roughly equal, consistent with a giant, but inconsistent with 
a supergiant. Add to that the lack of O I A7774 in J004618. 59+414410. 9, and we conclude 
these are also halo giants. 

We also include in Table the results of 6 other stars for which previous spectroscopy 
by one of us (N. C.) had identified strong O I A7774. Five of these are rank 1, and one is 
rank 2. (An additional two stars, one of rank 1 and one of rank 2, were in common with 
the WIYN spectroscopy, and agreed with our assessment.) These were all obtained with 
Hectospec using a lower resolution (5A) grating as part of a different project (Caldwell et 
al. 2009). In the following, we consider all stars found to have strong O I A7774 as certain 
supergiants, even if they were rank 2 based on their radial velocities, while we have removed 
the four stars shown spectroscopically not to be supergiants. 

It is of interest to see where our supergiants form in the various diagnostic diagrams we 
employed earlier. First, in Fig. [12]we superimpose the actual yellow supergiants on the CMD 
of M31. Note that indeed our bluest "yellow" supergiants extend into the blue supergiant 
region, as expected, given our relaxation of our original color selection to include stars as 
blue as B-V = (§[230. 

We next show in Figure [13] the locations of our yellow supergiants in the two-color 
diagram. We see that, had we relied upon a two-color plot to eliminate yellow supergiant 
candidates, we would have missed a number whose U — B colors are more negative than would 
be expected from the nominal supergiant sequence (shown in red). Of course, it could be 
that some of the stars with B — V > 0.4 but U — B < 0.1 will turn out to be foreground upon 
additional spectroscopy as none of the four rank 1 stars in this region yet have confirming 
O I A7774 spectroscopy. We have examined the spatial location of these four stars, however, 
and find that they are both in the spiral arms and near to other rank 1 yellow supergiants. 
we, therefore, feel that it is more likely this is a combination of photometric errors, slightly 
variable reddening, and uncertainties in the two-color relation for yellow supergiants. 

Finally, we show in Figure [H] the spatial distribution of our supergiant candidates 
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in M31, where blue symbols indicate the certain ones (either rank 1 or spectroscopically 
confirmed), and the red symbols the rank 2 candidates. For the most part, the yellow 
supergiants are found along the CO ring where star formation is most active, as expected. 

4. The H-R Diagram 

In Figure [15] we show the location of our 116 yellow supergiant candidates in the H-R 
diagram, along with the z = 0.040 evolutionary tracks. For simplicity, we have shown only 
the newer models which include an initial rotation of 300 km -1 ("S3" in Tabled]) as we view 
these as the most physically realistic. Several of these "S3" tracks were computed specifically 
for this project, and that of Massey et al. (2009). 

First, we find that the tracks do a good job of predicting the locations of yellow super- 
giants in the HRD. The most luminous yellow supergiants in our sample have log L/L & ~ 5.6. 
We do not find yellow supergiants with luminosities of (say) logL/L ~ 6, and this is in 
accord with what the evolutionary tracks predict. Note that with the older tracks (dotted 
tracks in Figure [3]) we might have expected to see some higher mass, warmer yellow super- 
giants. Indeed, the number of high luminosity yellow supergiants should have been similar 
to that seen for lower luminosity 12-2OM yellow supergiants, as the older 6OM track ex- 
tended into the yellow supergiant realm, and the duration of the yellow supergiant phase 
was 5,400 years (comparable to that of the lower masses as seen in Table [T]). 

This agreement with the new tracks is similar to what Massey et al. (2009) found for the 
coolest supergiants (T e g < 3800 K) in terms of the excellent agreement between the location 
of the tracks (and in particular the upper luminosities) and the observed locations of the 
stars. M31's RSGs have a maximum luminosity of log L/L Q ~ 5.4, a little bit lower than 
the most luminous yellow supergiants, as might be expected from the evolutionary tracks 
shown in Fig. [T5l as the 25M track does not extend to such cool effective temperatures. 

We now tackle the test we described in § 11.11 namely to see whether the relative number 
of yellow supergiants increases as we go to higher luminosities as the lifetimes in Table [1] 
suggest to us. A visual inspection of Figure [T5l says that the answer is clearly no: the number 
of yellow supergiants decrease. But, let us attempt a more quantitative assessment. 

We list in Table the number of yellow supergiants we find in each mass bin, both for 
the entire sample (rank 1 and rank 2) and for just the ones we are most certain are actually 
supergiants (rank 1). We then normalize these to the number of stars in the lowest mass 
bin, 12-15M Q . What we observe is a decreasing number of yellow supergiants as we go up 
in mass. At the bottom of Table E] we have included an "extra" mass bin, 15-25M , as it is 
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clear from Table [T] that the evolutionary tracks predict a much longer lifetime for the 2OM 
model than for either the 15M or the 25M , and we wanted to see what the agreement 
would look like if we ignored this track. 

We can estimate the number of yellow supergiants we expect from the models if we 
assume "steady state" star formation in M31. By this, we require only that the star formation 
rate averaged over the entire disk has stayed about the same for the past 20 Myr. In that 
case, the number of stars N in a particular evolutionary phase within a mass bin extending 
from one mass (mi) to another (m 2 ) will just be: 

NZ = [m r )Z x f 

where T is the slope of the initial mass function (taken here to be —1.35 following Salpeter 
1955; see also Massey 1998b), and f is the average duration of the evolutionary phase for 
masses mi and m 2 . In the final column of Table [6] we give the expected number according 
to the "S3" models. We have normalized the expected numbers relative to that of 12-15M . 

Here we find relatively poor agreement. According to the models, the evolutionary time 
scale for the yellow supergiant phase increases significantly with mass (Tabled]), more than 
compensating for the loss of stars due to the mass function. Thus we expect to find ~ 9 x 
more yellow supergiants between the 15-2OM tracks than between the 12-15M , or ~ 6x 
more between the 2O-25M tracks than between the 12-15M tracks. But, in reality, we find, 
0.7-0.8 x and 0.2 x as many, respectively. (Note that our results are insensitive to whether 
we count "all" of our candidates or just the certain ones.) Mostly this comes about because 
of the very long time predicted for the 2OM yellow supergiant stage relative to that of the 
lower masses (78,300 years vs 5,300 years). If we ignore the 2OM track we find only slightly 
better agreement, as the number of stars observed in the 15-25M is about the same as in 
the 12-15M track, while the models predict 4x as many. In addition, the models predict 
6x as many 25-4OM yellow supergiants as those of 12-15M . Based on this we expect to 
find 110-150 yellow supergiants with masses of 25-4OM , but instead we observe none. We 
can tell from Table [6] that the problem would be even worse if we had used the predictions 
of the "Old" z = 0.040 tracks, as the lifetimes are even longer for stars with masses >25M . 

We do note that the 4OM just barely enters the yellow supergiant realm. The long 
duration of this time stage reflects the fact that the star takes some time to turn around there. 
We also recall that in our sample we included stars that had bluer colors than our definition 
of yellow supergiants. So, if there were an abundance of such stars we would expect to have 
observed them, and, yet, none show up in the HRD. Nevertheless, if we were to assume that 
the duration of the 4OM yellow supergiant phase was years, rather than the 50,800 years 
we've included, then the number of expected stars between 25 and 4OM would be 1.5 x 
that of the 12-15M yellow supergiants. Thus rather than the 110-150 between the 25 and 



-21 - 



4OM we expect by using the 4OM lifetime, we would only expect 30-60 stars between 25 
and 40M Q . However, this still results in a significant discrepancy with observations, as we 
observe no yellow supergiants in this mass range. 

We emphasize that even though we expect some minimal contamination of our sample 
by foreground objects, the maximum contamination for the rank 1 objects ("mostly certain") 
is 15%. However, we obtain essentially the same ratios in Table [6] whether we count all of our 
candidates (rank 1 and rank 2) or just the rank 1 stars. So, it appears that our conclusion 
is robust. 

Recall also from Figure [8] that if anything our radial velocity survey was slightly biased 
towards the brighter stars than the fainter. We have made no allowance for this in the 
observed ratios in Table El but to do so would increase the discrepancy. The typical 12- 
15M stars have My = -7, or V = 17.8, while the 2O-25M stars have My = -9, or 
V = 15.8. The lower mass stars are therefore under represented by perhaps a factor of 
~ 1.5. Thus the 0.2 nominal observed ratio of the 2O-25M stars relative to the 12-15M 
stars should actually be lower by a factor of 1.5, suggesting that the disagreement with the 
5.6 ratio predicted by the models is about a factor of 40. 

Is this a problem with the higher mass evolutionary tracks predicting too long a time 
scale for the yellow supergiant stage, or with the lower mass tracks predicting too short a 
time scale? We can answer this indirectly by computing their expected lifetimes based upon 
the relative number of yellow supergiants and unevolved (O-type) stars observed. Using the 
LGGS data, Massey (2009) estimates the number of unevolved massive stars with masses 
> 2OM is about 24,800 in M31. The IMF-weighted H-burning lifetime is of order 5 Myr, 
and, assuming a constant star formation rate, we would thus expect to see 5 x 10~ 3 massive 
stars born each year. We observe 8 total (certain and probable) yellow supergiants above 
2OM . Recall that our sample contains only 68% of the stars located in the region for which 
we expect radial velocities to be < —150 km s -1 , and that region covered 73% of the area 
of the entire LGGS, from which the number of unevolved massive stars were estimated. We 
expect then that the true number of yellow supergiants with masses> 2OM is about 16. 
Therefore, we can estimate the actual ages of the yellow supergiant stage as 16/24800 x 5 
Myr. This is about 3200 years, which is consistent with the life times the 12-15M models 
predict, but is at variance with the much longer time scales predicted by the models for 
higher mass yellow supergiants. We suggest that these are too long by more than an order 
of magnitude. 

Could this discrepancy instead be an argument that the global star formation rate in 
M31 has in fact not been constant over the past 20 Myr? Yellow supergiants of 12M are 
roughly 17 Myr old, according to the models, while those of 25M are only 7 Myr. So, 
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if during that 10 Myr period the overall star formation rate had decreased by a factor of 
30-40 that would roughly compensate for the smaller number of stars that we find. However, 
such a drastic change is in conflict with other observations. Williams (2003) analyzed LGGS 
photometry and concluded that there has been a slight (25%) increase in the star formatation 
rate since 25 Myr ago. This is a marginal result, and consistent with constant star formation 
to within 2a (see his Table 2), but it certainly precludes the possibility of a 3000% - 4000% 
decrease over a similar time span. 



5. Summary, Discussion, and and Future Work 

We measured radial velocities for ~ 2900 stars in M31, identifying 54 as rank 1 (nearly 
certain supergiants) and 66 as rank 2 (probable supergiants). Follow up spectroscopy elim- 
inated 4 of the rank 2 stars, while confirming others as supergiants. The magnitude limits 
we chose should make our sample complete down to 12M . In all we observed 68% of the 
target candidates and the sample was restricted to the 73% of the LGGS area that should 
have radial velocities < —150 km s _1 . So, the true number of yellow supergiants should be 
about a factor of 2.0 larger than what we find. The foreground contamination proved to be 
96-98%. There may be a few halo yellow giants among our candidates, but comparison with 
the Besancon model suggests this should be minor, at most 15% for the rank 1 objects based 
purely upon the distribution of radial velocities. In practice we expect this contamination to 
be considerably less, since we replied upon the difference between the observed and expected 
velocities (where the latter is dependent upon position in the galaxy) to assign rank and 
membership. Nevertheless, it would be very useful to conduct follow up spectroscopy of the 
rank 2 objects in order to ascertain which have strong O I A7774 absorption or other spectral 
indicators of high luminosity. 

We compared the location and numbers of yellow supergiants in the H-R diagram to 
those expected from the Geneva evolutionary tracks. We find excellent agreement between 
the locations of stars in the H-R diagram and the tracks: there are not (for instance) high 
luminosity yellow supergiants with moderate temperatures that are inconsistent with the 
tracks. Rather, the inconsistencies we do note are related to the lifetimes predicted by 
the models for the yellow supergiant stage. The number of yellow supergiants decrease with 
increasing luminosity (mass), with no stars found more massive 25M . Yet, the long duration 
of the yellow supergiant phase predicted by the models for 20-40M Q suggests that we should 
see far more high luminosity yellow supergiants than what we observe. Comparing the 
number of yellow supergiants we find to the number of unevolved O-type stars, suggests that 
the typical duration of the yellow supergiant stage for stars with masses > 2OM should be 
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3000 years, similar to what the models predict for 12-15M Yet, the models predict lifetimes 
far greater than this. 

We do not yet have an adequate explanation for the discrepancy, but will address this 
in future work. The higher mass tracks (for which the predicted yellow supergiant lifetimes 
appear to be too long) show the stars evolving back to the blue after the RSG stage. If instead 
the stars ended their lives as RSGs, without this loop back to the blue, then the predicted 
yellow supergiant phase would be shorter as the star would pass through this region only 
once. This might be the case if the mass loss during the RSG stage had been significantly 
over-estimated. In part, this could be tested by comparing the number of observed Wolf- 
Rayet stars (WRs) with that of RSGs. The number of WRs in M31 is not known well enough 
to make this comparison as yet. Alternatively, it could be that the blue loops are present, 
but that the mass-loss rates during the evolution blue-wards have been underestimated. A 
sensitive test would be to conduct abundance studies of the yellow supergiants in this region 
of the HRD to look for evidence that any of these stars are in a post-RSG phase. An 
additional test would be to look for circumstellar material left from the slow dense wind of 
the RSG stage around any of these objects. 

It would, of course, be of interest to extend this work to galaxies with other metallicities, 
such as the Magellanic Clouds, where the unevolved massive star content is also known, 
and where the models predict long lifetimes for yellow supergiants even for 12-15M Q stars 
(Tabled]). We hope to carry out such work ourselves during the next observing season. 
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Fig. 1. — Color- magnitude diagrams. Upper: The LGGS photometry is used to construct a 
color-magnitude diagram for M31. Lower: The regions of large foreground contamination in 
the M31 diagram can be inferred from a CMD constructed using the Besancon model (Robin 
et al. 2003) of the Milky Way. 



- 28 - 



1 

Foreground 








: (aid 


% 




- 






















It 




_ , , , , i 




, ' , , i 





1 

B-V 




Foreground 
(Thick Disk) 



Sir 



1 

B-V 




Fig. 2. — Components of the foreground contamination of the M31 field identifed using the 
Besancon model. In the upper left we show the expected foreground contamination of the 
M31 field, where we have color-coded the giants (red), sub-giants (green), and the white 
dwarfs (cyan). The remaining points are primarily main-sequence dwarfs. In the other three 
panels we further break down the distribution into the disk component (upper right), the 
thick disk component (lower left) and the halo component (lower right). 
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Fig. 3. — The Geneva evolutionary tracks. The solid curves denote the latest models that 
include an initial rotation of 300 km s -1 , the dashed curves are the latest models that include 
no initial rotation, and the dotted curves are the older, non-rotating models. The various 
versions of the tracks appear in the same color for a given mass to reduce the confusion. The 
two vertical lines denote the yellow supergiant region, taken to be when the models have 
4800 < T c ff < 7500. The tracks are shown for three metallicities: z = 0.020 is characteristic of 
the solar neighborhood, z = 0.004 is characteristic of the SMC, and z = 0.040 is characteristic 
of M31. 
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Fig. 4. — Two-color diagram. The intrinsic two-color relationships are shown for dwarfs 
(green) and supergiants (red), from FitzGerald (1970). Upper: The points show the pho- 
tometry from the Local Group Galaxy survey of M31 (Massey et al. 2006). The supergiant 
sequence has been reddened for a typical E(B — V) = 0.13, with the reddening vector shown 
by the short line in the lower left. Lower: The points show the approximate foreground con- 
tamination from the Milky Way for the same solid angle and Galactic latitude and longitude 
as the M31 photometry. The data come from a simulation with the Besancon model (Robin 
et al. 2003). 
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Fig. 5. — The distribution of our sample across the face of M31. The M31 image is a mosaic 
of the 10 LGGS 36'x36' fields. The stars selected for our sample are shown by the red 
points. The requirement that the expected rotational velocity be < —150 km s -1 results 
in the "alligator jaw" pattern. The figure is slightly smaller than the area over which the 
LGGS has photometry, and hence a few points fall outside the image of the galaxy. 
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Fig. 6. — Examples of spectra and cross-correlations. We show sample spectra and the 
corresponding cross-correlation functions and fits. The thick bar shown on the spectra 
indicates the 4750-5550A, 5600-6500A regions used for the cross correlation. The upper 
star, J00303. 99+395939. 5, proves to be a foreground dwarf (§ 13.21) . with an average radial 
velocity of —48 km s" 1 compared to -534 km s -1 expected for its position. The lower star, 
J004156. 90+412109.0, proves to be an M31 supergiant (§ 13. 2p . with an average radial velocity 
of —335 km s _1 compared to —303 km s _1 expected for its position. Additional spectroscopy, 
discussed in § 13.41 confirms that the star has very strong 01 A7774, characteristic of a yellow 
supergiant. 
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Radial Velocities (km/s) from 17 Oct 2007 



Fig. 7. — Radial velocity constancy. We show here a comparison of the observed radial 
velocities for the two observations of the Brt5-1 field. The line shows the 1:1 relation. 
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Fig. 8. — Distribution in magnitude. We compare the magnitudes of the subsample of 
stars for which we obtained radial velocities (red, dashed histogram) with that of the parent 
population (black, solid histogram). 
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Fig. 9. — Comparison between the observed and expected radial velocities. Here we plot 
the difference between the observed and expected radial velocities vs. the expected radial 
velocities. Any M31 members will lie near a difference of 0, while the dark band is composed 
of foreground dwarfs. Those objects with V Q b s < -400 km s _1 are marked in green and those 
with -400 km s~ x < V Q b s < -300 km s _1 as red. The assigned "ranks" (1= mostly certain 
supergiant; 2=probable supergiant) are shown. 
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Fig. 10. — Histogram of radial velocities. The observed radial velocities (red, dashed his- 
togram) is compared with what we expect for the foreground contamination (black, solid 
histogram) computed from the Bescancon model. 




Fig. 11. — Comparison of effective temperature scales. The Atlas9 model predictions are 
shown by open (black) circles for surface gravities appropriate to supergiants. The filled (red) 
circles are from the Kovtyukh (2007) study of F- and G-type supergiants. The solid (green) 
curve is the Flower (1996) relation for supergiants. the two horizontal lines correspond to 
effective temperatures of 7500 K and 4800 K, the region we are considering the realm of the 
yellow supergiants here. 
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Fig. 12. — Color-magnitude diagram revisited. Same as Figure [TJ but the addition of the 
newly confirmed yellow supergiants marked as colored filled circles. 
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Fig. 13. — Two-color diagram revisited. Same as Figure H] left, but with the yellow super- 
giants from this paper now marked by filled circles (if rank 1 or spectroscopically confirmed 
rank 2) or as x's (if unconfirmed rank 2). 
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Fig. 14. — The distribution of our yellow supergiants across the face of M31. The blue 
symbols represent the rank 1 (certain) yellow supergiant candidates, while red represents 
the rank 2 (less certain) candidates. Compare with Figure [51 
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Fig. 15. — H-R diagram. For simplicity, we show only the latest (z = 0.040) Geneva tracks 
with an initial rotation of 300 km s _1 ; the location of the older tracks can be seen in Fig- 
ure [31 The solid points are our certain yellow supergiants (either rank 1 or spectroscopically 
confirmed rank 2) while the open points are the as- yet unconfirmed rank 2 yellow supergiant 
candidates. 
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Table 1. Theoretical Yellow Supergiant Duration (yr) 





Solar Neighb 
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Mass 




z = 0.020 




z 
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= 0.040 


(M ) 


S3 b 


so b 


01d c 


S3 d 
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56500 
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5400 
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403600 


198100 


148100 


32700 


188200 


50800 
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25 


5400 


73000 


23600 
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18500 




666000 


20 
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92500 


64800 
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16500 
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3800 


15 
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206600 


60400 


92000 


3200 




1900 


12 


3600 


32800 


51000 


33200 


30700 


21800 


5300 
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a For the purposes of this calculation, yellow supergiants are defined as having effective 
temperatures between 4800 K and 7500 K. Note that not all versions of the evolutionary 
models are available for each mass and metallicity. 

b S3 has initial rotation 300 km s -1 , and SO has an initial rotation of km s -1 ; ages were 
determined using models from Meynet & Maeder (2003). 

c Rotation not included; determined using models from Schaller et al. (1992). 

d S3 has initial rotation 300 km s" 1 , and SO has an initial rotation of km s _1 ; ages were 
determined using models from Meynet & Maeder (2005), except for the 9, 12, 20 and 25 
Mq, z = 0.040 models newly computed for this study and that of Massey et al. (2009), 
and the 12, 15 and 20 M , z = 0.004 models from Maeder & Meynet (2001). 

c The SO models for z = 0.004 were computed with an initial rotation of km s -1 ; ages 
were determined using the models from Maeder & Meynet (2001). 

f Rotation not included; determined using models from Charbonnel et al. (1993). 

g Rotation not included; determined using models from Schaerer et al. (1993). 



Table 2. Configurations Observed 



Config 0:2000 £2000 Exps UT Date(s) N stars observed 
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Table 3. Stars with Observed Radial Velocities* 



Star 


a 2000 


•52000 


Vcl oba 


r a 


Vel expect 


Vcl oba -Vcl oxp 


V 


B - V 


Rank b 








km s~ 1 




km s — 1 


km s — 1 








J003702. 13+400945. 6 


00:37:02.127 


+40:09:45.53 


14.8 


19.8 


-513.5 


528.3 


16.09 


1.29 


3 


J003702. 47+401742. 5 


00:37:02.467 


+40:17:42.42 


-26.4 


21.3 
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1.32 


3 


J003702. 80+400516. 8 
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-4.5 
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0.71 
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32.6 
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0.84 
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3 



*Thc full version of this tabic is available in the on-line edition. 
a Tonry Ik. Davis 1979 r parameter. 

b R.ank: 1 — highly likely supergiant; 2— probable supergiant; 3— dwarf; CI— cluster 
c Radial velocity observed on two different nights differed by 10 km s — 1 or more. 
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Table 4. Observed Properties of Probable M31 Members 
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3 


19. 


2 


-163 


5 


-24 


8 


15. 


60 





.46 


2 


.7004409.23+415941 


1 


00: 


44 


:09 


227 


+41 


:59: 


40. 


93 


-102 





37. 


1 


-156 


1 


54. 


1 


17 


42 





.62 


2 


.7004409.98+420121 


1 


00: 


:44 


:09 


977 


+42 


:01: 


20. 


92 


-145 


1 


23 


5 


-159 





13. 


9 


17 


42 





.60 


2 


.7004410.62+411759 


7 


00: 


44 


10 


617 


+41 


:17: 


59. 


57 


-244. 


5 


21. 


1 


-243 


3 


-1. 


2 


17. 


63 





.40 


2 


Mag-253496 




00: 


:44 


12 


450 


+41 


:16: 


08. 


42 


-226 


9 


54. 


5 


-253. 


6 


26. 


7 


15 


15 


1 


.18 


2 


.7004424.21+412116. 





00: 


:44 


24 


207 


+41 


:21: 


15. 


86 


-238. 





14. 


3 


-225 


9 


-12. 


1 


16. 


73 





.91 


2 


.7004427.76+412209 


8 


00: 


44 


:27 


757 


+41 


:22: 


09. 


66 


-213 


9 


19. 


9 


-221 


3 


7. 


4 


17. 


26 


1 


.03 


2 


.7004428.99+412010 


7 


00: 


:44 


.28 


987 


+41 


:20: 


10. 


56 


-216 


8 


9. 


1 


-233 


9 


17. 


1 


17. 


80 





.31 


2 


.7004432.01+412442 





00: 


44 


32 


007 


+41 


:24: 


41. 


86 


-168 


2 


5. 


4 


-205. 


3 


37 


1 


18. 


28 





.07 


2 


.7004432.41+412947 


5 


00: 


:44 


32 


407 


+41 


:29: 


47. 


35 


-120 


2 


59. 


4 


-159 


3 


39 


1 


17. 


05 





.73 


2 


.7004440.60+412704 


1 


00: 


44 


40 


597 


+41 


:27: 


03. 


96 


-189 


8 


27 


8 


-192 


7 


2. 


9 


18. 


22 


1 


.33 


2 


J004441. 56+412636 


6 


00 


:44 


41 


557 


+41 


:26: 


36. 


46 


-165 


4 


9 


3 


-196 


9 


31. 




17. 


96 





.44 


2 


J004444. 50+412314 


3 


00 


:44 


44 


497 


+41 


:23: 


14. 


16 


-186 


3 


40. 


6 


-220 


8 


34. 


5 


17. 


57 





.65 


2 


J004447. 45+412409. 


7 


00 


:44 


47 


447 


+41 


:24. 


09. 


56 


-197 


9 


12. 


3 


-216 


4 


18. 


5 


18. 


10 





.22 


2 


J004458. 01+413217 


5 


00 


:44 


:58 


007 


+41 


:32. 


17. 


35 


-128 


8 


11. 


6 


-165 


7 


36. 


9 


17. 


79 





.17 


2 


J004508. 90+413117 


8 


00 


:45 


08 


897 


+41 


31 


17. 


65 


-144. 





32. 


6 


-182 


5 


38 


5 


17. 


12 





.61 


2 


J004509. 86+413031. 


5 


00 


:45 


09 


857 


+41 


:30: 


31. 


35 


-402 


6 


6. 


3 


-188 


3 


-214. 


3 


16. 


83 





.56 


2 


J004518. 17+413615. 


6 


00 


:45 


18 


167 


+41 


:36: 


15. 


45 


-101 


1 


41. 


2 


-154 


4 


53 


3 


17. 


84 





.68 


2 


J004518. 76+413630. 


7 


00: 


45 


:18 


757 


+41 


:36: 


30 


55 


-128 


4 


19. 


5 


-153 


1 


24. 


7 


16. 


70 





.51 


2 


J004526. 93+412613. 


6 


00: 


45 


26 


927 


+41 


:26: 


13. 


46 


-255 





24. 


6 


-218 


8 


-36 


2 


18. 


45 





.75 


2 


J004532. 62+413227. 


8 


00 


45 


:32 


617 


+41 


:32: 


27. 


65 


-436. 


7 


35. 


3 


-190 


2 


-246. 


5 


15. 


79 





.85 


2 


J004535. 23+413600. 


5 


00 


45 


:35 


227 


+41 


:36: 


00 


35 


-154 


3 


75. 


8 


-170 


8 


16 


5 


15. 


88 





.94 


2 


J004554. 48+413359. 


8 


00: 


45 


:54 


.477 


+41 


:33: 


59. 


65 


-344 


4 


20. 


1 


-193 


1 


-151. 


3 


17. 


16 





.59 


2 


J004559. 84+414038. 


2 


00: 


45 


:59 


837 


+41 


'40: 


38. 


04 


-161 


9 


5. 


9 


-161 





-0 


9 


17 


71 





.22 


2 


J004618. 59+414410 


9 


00: 


:46 


:18 


587 


+41 


44: 


10. 


74 


-149 


2 


34. 


8 


-154 


6 


5. 


4 


15. 


22 





.65 


2 


J004658. 64+414948. 


4 


00: 


46 


:58 


637 


+41 


:49: 


48. 


24 


-134 


6 


15. 





-153 


1 


18. 


5 


18 


30 





.48 


2 


Mag-237751 




00 


:41 


01 


181 


+41 


13: 


45. 


83 


-441. 


5 


41. 


1 


-348 





-93 


5 


14. 


13 


1 


.12 


CI 


J004345. 23+410608. 


5 


00 


43 


45 


227 


+41 


06: 


08. 


38 


-402 


1 


28. 


3 


-298 


2 


-103. 


9 


18. 


32 





.52 


CI 


J004356. 46+412203. 


3 


00 


43 


56 


.457 


+41 


:22: 


03. 


16 


-362 





59. 


1 


-203 


1 


-158. 


9 


16. 


88 





.08 


C7 


J004358. 15+412438. 


8 


00 


43 


:58 


147 


+41 


24: 


38 


66 


-380. 


4 


47. 


9 


-171 


4 


-209. 





17. 


35 





.99 


CI 


,7004403.98+412618. 


7 


00 


:44 


03 


977 


+41 


26: 


18 


56 


-293 


1 


26. 


5 


-157 


6 


-135. 


5 


18. 


17 





.71 


C7 


J004446. 42+412918. 


3 


00 


:44 


46 


417 


+41 


29: 


18 


16 


-209 


8 


63. 


5 


-179 


4 


-30 


4 


17. 


96 





.26 


C7 


J004545. 58+413942. 


4 


00 


45 


45 


.577 


+41 


:39: 


42. 


25 


-204 


9 


48. 


5 


-155 


5 


-49 


4 


18. 


13 





.89 


C7 



a Tonry Sz Davis 1979 r parameter. 

b Rank: 1 — highly likely supergiant; 2— probable supergiant; 3— dwarf; CI— cluster 
c Radial velocity observed on two different nights differed by 10 km s — 1 or more. 



-46- 



Table 5. Derived Properties of Potential M31 Supergiants 





Star 




Rank 


My 


T .ff 


log (L/Z, ( 


Mag-253496 




2 


(-9.65) 




71 ^"i 




° ' J 


J004251 


.90+413745 


.9 


2 


(-9.82) 




792) 


y° 


84) 


J004618 


.59+414410 


.9 


2 


(-9.58) 




787) 


(.*>• 




J004406 


.32+420131 


.3 


2 


-9.20 


3 


822 


5. 


■58 


J004532 


.62+413227 


.8 


2 


-9.01 


3 


758 


5. 


55 


J004535 


.23+413600 


.5 


2 


-8.92 


3 


746 


5. 


53 


J004034 


.82+401825 


.5 


2 


-8.50 


3 


794 


5. 


31 


J004247 


.30+414451 


.0 


2 


-8.39 


3 


821 


5. 


25 


J004318 


.57+415311 


1 


2 


-8.26 


3 


768 


5. 


24 


J004424 


.21+412116 





2 


-8.07 


3 


.750 


5. 


18 


J003907 


.59+402628 


.4 


1 


-8.06 


3 


757 


5. 


17 


J004021 


.21+403117 


.1 


1 


-8.15 


3 


867 


5. 


15 


J004518 


.76+413630 


.7 


2 


-8.10 


3 


812 




14 


J004129 


.31+405102 


.9 


1 


-7.93 


3 


776 


5. 


10 


J004509 


.86+413031 


.5 


2 


(-7.97) 


( 3 


804) 


\° 


09) 


J004207 


.22+405148 


.3 


1 


-7.77 


3 


743 


5. 


08 


J004337 


.16+412151 


.0 


2 


-7.78 


3 


750 


5. 


07 


J004207 


.85+405152 


.4 


1 


-7.81 


3 


771 


5. 


06 


J003930 


79+401841 


1 


1 


-7.90 


3 


823 


5. 


05 


J004313 


.02+414144 


.9 


2 


-7.83 


3 


791 


5. 


.04 


J004101 


.24+410434 


.6 


1 


-7.68 


3 


749 


5. 


03 


J004432 


.41+412947 


.5 


2 


-7.75 


3 


776 




03 


J004120 


.56+403515 


.4 


1 


-7.70 


3 


796 


4 


99 


J004427 


.76+412209 


.8 


2 


-7.54 


3 


734 


4 


99 


J004508 


.90+413117 


.8 


2 


-7.68 


3 


794 


4 


98 


J004248 


.11+403434 


.5 


2 


-7.59 


3 


771 


4 


97 


J004243 


.54+414620 


.6 


2 


-7.60 


3 


774 


4 


97 


J003949 


86+405305 


.8 


1 


-7.46 


3 


736 


4 


96 


J004554 


.48+413359 


.8 


2 


-7.64 


3 


797 


4 


96 


J004156 


.90+412109 


.0 


1 


-7.63 


3 


807 


4 


96 


J003745 


.26+395823 


.6 


1 


-7.64 


3 


797 


4 


96 


J003943 


.43+403524 


.9 


1 


-7.63 


3 


879 


4 


95 


J003953 


.55+402827 


.7 


1 


-7.35 


3 


712 


4 


95 


J003934 


.02+404714 


.2 


2 


-7.59 


3 


795 


4 


94 


J004259 


.95+410220 


.3 


1 


-7.39 


3 


730 


4 


94 


J004101 


.76+410429 


.2 


1 


-7.38 


3 


762 


4 


89 


J004237 


.52+413024 


.8 


2 


-7.47 


3 


877 


4 


89 


J004232 


.90+412103 


.4 


2 


-7.37 


3 


756 


4 


89 


J003948 


.85+403844 


.8 


1 


-7.40 


3 


787 


4 


88 


J004409 


.98+420121 


1 


2 


-7.38 


3 


796 


4 


86 


J004409 


.23+415941 


1 


2 


-7.38 


3 


793 


4 


80 


J004038 


.10+403827 


.2 


1 


-7.21 


3 


743 


4 


85 


J003937 


.44+394941 


1 


2 


-7.30 


3 


823 


4 


82 


J004107 


.40+405328 


.6 


2 


-7.30 


3 


821 


4 


81 


J004444 


.50+412314 


.3 


2 


-7.23 


3 


787 


4 


81 


J004009 


.13+403142 


.9 


1 


-7.17 


3 


760 


4. 


81 


J004410 


.62+411759 


.7 


2 


-7.17 


3 


835 


4. 


76 


J004229 


.30+405727 


.6 


1 


-7.17 


3 


830 


4. 


76 


J004301 


.96+405315 


.2 


2 


-7.02 


3 


.767 


4. 


74 


J004559 


.84+414038 


.2 


2 


-7.09 


3 


884 


4. 


74 


J004217 


.15+403740 


.2 


2 


-7.03 


3 


.787 


4. 


73 


J004255 


.16+413515 


.3 


2 


-6.96 


3 


761 


4. 


73 


J004034 


.00+403500 


.1 


1 


-7.10 


3 


.827 


4. 


73 


J003935 


.23+401947 


.7 


1 


-6.88 


3 


741 


4. 


72 


J003948 


.28+403856 


.4 


1 


-6.92 


3 


.750 


4. 


72 


J004458 


.01+413217 


.5 


2 


-7.01 


3 


898 


4. 


71 


J004440 


.60+412704 


.1 


2 


-6.58 


3 


688 


4. 


70 


J004250 


.67+414008 


.6 


2 


-6.95 


3 


785 


4. 


70 


J004302 


.55+413332 


.0 


2 


-6.93 


3 


.773 


4. 


70 


J004518 


.17+413615 


.6 


2 


-6.95 


3 


.783 


4. 


70 


J004128 


,74+405224 


.7 


1 


-7.01 


3 


876 


4. 


70 


J004428 


.99+412010 


.7 


2 


-7.00 


3 


857 


4. 


69 


J004021 


.64+403256 


.5 


1 


-6.97 


3 


813 


4. 


69 


J004235 


.91+413050 


.6 


2 


-7.00 


3 


820 


4. 


69 


J003926 


.72+402239 


.4 


1 


-6.73 


3 


961 


4. 


67 


J004149 


.87+412712 


.7 


1 


-6.81 


3 


.761 


4. 


66 


J004025 


.48+405041 


.0 


1 


-6.81 


3 


932 


4. 


66 



Comment 



No O I A7774-Nonmcmbcr 
No O I A7774-Nonmember 
Non- member 



Strong O I A7774 
Very strong O I A7774 
Strong O I A7774 

Weak O I A7774-Non member? 
Strong O I A7774 



Very strong O I A7774 



Very strong O I A7774 
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Table 5 — Continued 





Star 




Rank 


M v 


T, 


= // 


log(£ 


/L 


J004214 


85+405652. 





1 


-6.83 


3 


.891 


4. 


G4 


J004131 


50+403917. 


8 


1 


-6.84 


3 


803 


4. 


G4 


J004304 


.89+410345 


9 


1 


-6.73 


3 


753 


4. 


G4 


J004441 


.56+412636 


6 


2 


-6.84 


3 


827 


4. 


63 


J004245 


.97+414422. 


1 


2 


-6.81 


3 


896 


4. 


63 


J004314 


47+414229 


1 


2 


-6.67 


3 


951 


4. 


63 


J004032 


90+404352. 


8 


1 


-6.68 


3 


749 


4. 


63 


J004233 


76+410014 


6 


1 


-6.77 


3 


783 


4. 


63 


J004002 


91+400659. 


2 


2 


-6.82 


3 


836 


4. 


62 


J004240 


34+412807. 


2 


2 


-6.79 


3 


810 


4. 


G2 


J004034 


00+405358. 


3 


1 


-6.78 


3 


823 


4. 


61 


J004143 


45+403956. 


4 


1 


-6.80 


3 


833 


4. 


61 


J003711 


98+395445. 


2 


2 


-6.74 


3 


803 


4. 


60 


J003930 


55+403135. 


2 


1 


-6.65 


3 


769 


4. 


59 


J004432 


01+412442 





2 


-6.52 


3 


966 


4. 


59 


J004311 


.34+414240 


9 


2 


-6.63 


3 


767 


4. 


59 


J004007 


14+410321 


8 


2 


-6.72 


3 


807 


4. 


59 


J004035 


37+405701 





1 


-6.69 


3 


888 


4. 


58 


J003936 


96+400743. 


8 


2 


-6.61 


3 


765 


4. 


58 


J004144 


76+402808. 


9 


1 


-6.51 


3 


734 


4. 


58 


J004118 


69+403152 





1 


-6.66 


3 


.784 


4. 


58 


J004447 


45+412409 


7 


2 


-6.70 


3 


882 


4. 


58 


J004247 


25+410039 


2 


1 


-6.59 


3 


765 


4. 


57 


J004252 


87+412328. 


6 


2 


-6.50 


3 


954 


4. 


57 


J004029 


38+403604. 


2 


1 


-6.67 


3 


.857 


4. 


56 


J003942 


35+404031. 


8 


2 


-6.65 


3 


806 


4. 


56 


J004120 


99+403453. 


5 


1 


-6.67 


3 


825 


4. 


56 


J004338 


76+414915 


1 


2 


-6.64 


3 


871 


4. 


55 


J004303 


69+414543 


3 


2 


-6.59 


3 


884 


4. 


54 


J004030 


10+402943. 


1 


1 


-6.53 


3 


924 


4. 


54 


J004053 


97+403256. 


2 


1 


-6.49 


3 


766 


4. 


53 


J004325 


93+413910 


6 


2 


-6.48 


3 


773 


4. 


52 


J004212 


20+405513. 


9 


1 


-6.56 


3 


881 


4. 


52 


J004020 


37+410723. 


2 


2 


-6.48 


3 


783 


4. 


51 


J004348 


01+415406. 


2 


2 


-6.50 


3 


797 


4. 


51 


J004009 


70+403719. 





1 


-6.53 


3 


889 


4. 


51 


J004658 


64+414948 


4 


2 


-6.50 


3 


819 


4. 


50 


J004223 


21+412803. 


2 


2 


-6.43 


3 


785 


4. 


49 


J004030 


62+404523. 


8 


1 


-6.47 


3 


864 


4. 


48 


J004017 


72+400436. 


6 


1 


-6.47 


3 


873 


4. 


48 


J004102 


78+410900 


6 


1 


-6.45 


3 


806 


4. 


48 


J004215 


06+405148 


3 


1 


-6.32 


3 


.754 


4. 


48 


J003922 


08+402031. 


5 


1 


-6.44 


3 


897 


4. 


48 


J004226 


53+410123. 


9 


1 


-6.43 


3 


820 


4. 


47 


J004526 


93+412613 


6 


2 


-6.35 


3 


.772 


4. 


47 


J004231 


84+412039. 


9 


2 


-6.36 


3 


804 


4. 


45 


J003955 


87+401636. 


4 


2 


-6.37 


3 


807 


4. 


45 


J004026 


68+403604. 


4 


1 


-6.36 


3 


.878 


4. 


44 


J004020 


37+403428. 


8 


1 


-6.37 


3 


876 


4. 


44 


J004148 


69+413814 


1 


2 


-6.31 


3 


790 


4. 


44 


J004201 


09+403951. 


9 


1 


-6.34 


3 


.845 


4. 


43 


J004259 


.62+413946. 


1 


2 


-6.32 


3 


816 


4. 


43 


J003725 


57+400731. 


9 


2 


-6.33 


3 


815 


4. 


43 



Strong O I A7774 
Very strong O I A7774 
Strong O I A 7774 



Strong O I A7774 



Strong O I A7774 
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Table 6. Number of Yellow Supergiants Compared to Models 



Mass 
Range 


# 
All 


# 

Certain 


Ratio relative to 12-15M 


All 


Certain 


Models 


12-15M© 


41 


20 


1.0 


1.0 


1.0 


15-2OM 


28 


16 


0.7 


0.8 


8.7 


20-25M© 


8 


4 


0.2 


0.2 


5.7 


25-40M© 








0.0 


0.0 


5.5 


15-25M© 


36 


20 


0.9 


1.0 


3.6 




Radial Velocities (km/s) from 17 Oct 2007 



